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I. GENERAL DISCUSSION 

The purpose of this article is to analyze the dependence of the energy of an 
elementary excitation on the strength of the confinement potential, which exists 
in a planar semiconductor heterostructure. Due to the fascinating technological 
progress in the field of man-made structures, it has become possible to fabricate e.g. 
quantum wells of a widely varying shape. It is an interesting theoretical task to 
. discuss the excitation spectrum of such semiconductor structures as function of the 

tunable parameters, such as well width, well height, etc.. Concerning the excitations 
of interest, we concentrate on particle-phonon systems, the particles being electrons 
^ ; or holes. The simplest example is that of a single polaron, that is an electron, 

coupled to a certain branch of lattice vibrations. Another example is that of a 
polaronic exciton, that is an electron-hole pair, coupled to phonons. Whereas the 
latter one is important to characterize optical properties, the former one has direct 
implications for the transport behavior of the materials of interest. 

We assume that the interface(s)-induced confinement can be mimicked by a sim- 
ple potential V n (z n ), n being the particle number, z n the corresponding coordinate 
(the growth direction of the heterostructure will always be assumed as z-direction) . 
Explicit forms of V n (z n ) may be rectangular wells, parabolas etc.. In addition, we 
suppose translation invariance to hold within the xy-plane. We remark that effects 
as surface roughness would destroy this property and could lead to the appearanc 
of new phenomena (e.g. localized states). 
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In the following equation, we define the class of models under discussion: 
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The nomenclature is self explaining. The quantity L r (r 1 ,r 2 ) is to contain the con- 
finement potentials as well as the particle interaction: 
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U(ri,..,r N ):=J2 V n{zn) + 2 H V n y{r n , r n /), (2) 

n, i = 1, 
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where V n ,n' has to be calculated as potential energy of particle n, exposed to the elec- 
trostatic potential of particle n' . Because of the boundary conditions, V n)n i(r n , r n i) 
itself is not translation invariant (see e.g. Ref. |1|]). The particle-phonon coupling is 
of Frohlich type. The most prominent example to be used here is that of a coupling 
to (LO)-phonons. 

The model has two relevant limiting cases, which should be reproduced by any 
theory. Let the maximum of the well widths be L and the minimum V . If V tends 
to infinity, the confinement is irrelevant and the energy spectrum of H is that of a 
three-dimensional well-material excitation. If L tends to zero, the (finite height) well 
is irrelevant, leaving us with the spectrum of a three-dimensional barrier-material 
excitation. The behaviour for intermediate values of the well widths can qualita- 
tively be discussed as follows. Varying L, V from sufficiently large values to smaller 
ones, the binding energy should increase due to the higher Coulomb correlation (for 
instance, the reader should recall that the energy of the two-dimensional hydrogen 
ground state is four times larger than that of a three-dimensional one). When L, V 
become smaller and smaller, the ground-state wave function will more and more 
effectively tunnel into the barrier material — the energy approaches the barrier 
limit. 

Thus, we might expect a maximum of the binding energy to appear at interme- 
diate values of L, L'. It was a controversely discussed question whether or not this 
maximum appears at relevant (that is not too small) values of L. The answer to 
this question might be not the same for different systems. 

II. POLARONS 

The physics of polarons, confined to quantum wells, passed a few stages, and 
it is not possible to present here even a brief list of references. In particular, it 
was found that different phonon modes contribute to the polaron binding energy - 
confined bulk 2phonons inside the well, interface phonon mode and half-space bulk 
phonon mode in the barrier. We cite only papers concerning polarons confined 
to a finite rectangular potential (one layer heterostructure) where contribution of all 
phonon modes were taken into account. Anyway, there are problems to be addressed 
while dealing with multilayered heterostructures. Namely, we have to answer the 
following questions: 

1) How to deal with multilayered heterostructures? The total number of phonon 
modes becomes too large to make numerical calculations even with modern com- 
puters. Besides, a multilayered heterostructure can generate a confining potential 
of rather complicated form, not only the rectangular one. 

2) How to deal with mass- and dielectric mismatches in different layers? The 
polaron effective mass m(z), the electron-phonon coupling constant a(z) and the 
phonon dispersion law do depend on a layer, that is, on the electron position. To 
glue solutions in different layers seems to be a cumbersome job. 
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To tackle these problems we suggest specific approximations, which will briefly 
be indicated here. 



• A multilayered GaAs / Al x Gai_ x As heterostructure is considered as an effec- 
tive medium. Its mean parameters are to be defined by averaging over different 
layers according to the way they enter the Hamiltonian. 

• The bulk phonon mode only inhabits an effective medium with mean charac- 
teristics. 

We specify the electronic part of the Hamiltonian: 

H el = H el>ll + H d>± = |L + it + V (z), (3) 

2m 2m 

The mean electron band mass m is defined by the equation 

flW,±^ = ^i, -= Jdz 1 -^^, (4) 

m J m(z) 

— oo 

where ipi(z), E% are the ground state wave function and the energy for the electron 
motion in z direction. As ipi and E\ depend on m, we actually have the system of 
two equations @ to calculate the mean band mass m. 
The free LO-phonon Hamiltonian reads as follows: 

oo 

Hph = fwj^o 4 a it ' Ulo= J dz u ( z ) \^( z )\ 2 ■ ( 5 ) 

k — oo 

As m is found already, we define here the mean phonon frequency uj lo . Note that in 
this paper we are not interested in processes of emission, absorption or scattering of 
phonons. Instead we concentrate on virtual phonons in a cloud around an electron. 
Subsequently, the properties of the effective phonons do depend on the position of 
the electron as it follows from Eq. (||). 

In the same way we define the effective electron-phonon interaction Hamiltonian 
in the standard Frohlich form with the mean Frohlich coupling constant a: 



2 u{z) ( , , / mw LO 
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a = / dz a(z) L " . (6) 

J uj lo V y m(z)uj{z) ) 

— OO ' 

Evidently, this model belongs to the class defined in Eq. (|l|) As examples we 
studied 1) a one-layer heterostructure described by a rectangular confining potential 
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(the z-dependence of the masses and dielectric parameters is completely analogous) 
and 2) a multilayered heterostructure corresponding to the Rosen-Morse potential 

V(z) = V tanh 2 (V— ) . (8) 

We use perturbation theory in powers of a for both potentials, but in the first 
case we perform the summation over all virtual states while in the case of the 
Rosen- Morse potential the Green function (see HH) can be used. To compare re- 
sults for the Rosen-Morse and the rectangular potentials, an effective width L of 
the Rosen-Morse potential has to be found. We define it as the width of a rect- 
angular potential of the same height Vo with the same ground-state energy. The 
dependence L(Lrm) can then be calculated. The parametrization for experimental 
data concerning GaAs / Al x Ga\_ x As heterostructure is based on the results reported 
in Ref. || with some modifications, which are discussed in our paper 0. Actually 
we use the dependence of the parameters on the Al mole fraction x which depends 
in turn on the coordinate z via the relation V(z) = 600 • (1.155a; + 0.37x 2 ) meV. 
The confining potential V(z) being given, we know the dependence x(z) and, sub- 
sequently, the values of the parameters a, m, uj at each point of the heterostructure 
which are averaged then following Eqs. (f|), (|5|) and 
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FIG. 1. The polaron binding energy and the effective mass in the rectangular poten- 
tial. Contributions of the discrete Adi s E (A^m/m) and continuous A con E (A con m/m) 
spectrum are shown as well as the so called leading term approximation A^E (Aum/m) 
when only the ground state is taken into account as an intermediate virtual state (dashed 
line). 

The polaron energy and effective mass are calculated for x = 0.3. Peaks are 
found for the effective mass at some potential widths, while the energy demonstrates 
rather a smooth behavior between the correct 3D-limits as is seen in Fig. |. As to 
the Rosen-Morse potential, the results are presented in Fig. |2| together with those 
for the rectangular potential of the corresponding effective width. One can see an 
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excellent coincidence of the results obtained within the different techniques; clearly, 
this fact increases their reliability. A comparison is also made with the results of 
the papers @|§, and the details are discussed in our paper QJ. 
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FIG. 2. Our results for the Rosen-Morse potential in comparison with these for the 
rectangular potential (dashed line) of the same effective width L{Lrm)- 



III. EXCITONS 

Sampling the previous literature, most work has been done on rectangular quan- 
tum wells with confinement potentials of type fl7|). The electron-hole potential can 
be calculated as indicated above and was given e.g. in Ref. [TJ. 

To treat eigenvalue problems as the present one, we use tractable decompositions 
of the Hamiltonian to generate lower bounds for the ground-state energy. The basic 
idea is as follows: Assume we study the Hamiltonian H = p 2 z /2m + Vi(z) + V 2 (z) to 
find its ground-state energy E. Then we use the decomposition 

2 2 

H 1 = x^ + V 1 {z) 1 H 2 = (l-x)^- + V 2 (z), 0<x<l. (9) 

If Ei(x),E 2 (x) are the corresponding ground-state energies of Hi,H 2 , then a lower 
bound for E is: E > max x (Ei(x) + E 2 (x)). 

Upper bounds are produced by variational methods: The trial wave-function 
used in our calculations had the form: 

*(*=!, *i,*2) = $i(^i)$2(^)e- a v / -i+^-^) 2 , (io) 

where $i(zi) are the ground-state eigenfunctions of the free electron (i = 1) or 
the hole (i = 2) in the confining potentials of the type (J]). Evidently, the varia- 
tional parameters a, b can be used to fit 3D and 2D limiting cases. If the masses 
can be assumed as constant over the heterostructure, these methods can profitably 
be combined with functional-integral techniques. Fig. |] shows our result || for 
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Alo.3Gao.7As/GaAs/Alo.3Gao.7As. in comparison with experimental |10|Jll"H and pre- 
vious theoretical results |§. Clearly, the maximum appears at a relevant width. 




FIG. 3. Comparison of results for the binding energy of an exciton in a rectangular 
quantum well as function of the well width. 

A second class of confinement potentials is of parabolic type, that is 

V5W = I") 

where A» denotes the dimensionless confinement strength, i?^ is the Rydberg unit, 
which was extracted for reasons of convenience. To study the confinement-induced 
effects on the spectrum as accurately as possible, we disregarded any parameter 
mismatch. The quantity of interest is the diagonal element of the reduced density 
operator, namely 

Pp{C) :=< C\tr Ph e-P H \C > . (12) 

In this formula C is an abbreviation for an arbitrary (but fixed) set of the posi- 
tion coordinates of the particles involved. The right-hand side of Eq. (|1^) can be 
represented functional integral 

Pp{C) = Z Ph J 5 6 R e' 3 ™. (13) 

In Eq. flTBD Zph is the free-phonon partition function, and S reads as follows: 



S[R] := f Q dr ^R^(r) + U^r), R 2 (r))) 



2 $ $ 

£ £ j j dTdT , G{j _ r /) e ik[R„(r)-R„,(r')]_ (14) 



n,n'=l k 
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Within the functional integral (|13D, / 5 6 R.... is to indicate integration over all real, 
6- dimensional paths R(t), which start and end at the point C. The kernel function 
G(t — t') is defined as 

_Rw(j3-|r|) _|_ p hw\r\ 

G(t) := n^-—-^ • (15) 

v > 2 [eP Jiu - 1] 

It is well known that functional integrals of type ( 13|) with an action (li~4[) cannot 



be evaluated in analytical form. Starting from the exact expression, we use varia- 
tional procedures as in Feynman's famous paper on polarons to find upper bounds 
on the ground-state energy. The necessary input is a trial action, which is accessible 
to a numerical treatment. 
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FIG. 4. Binding energy of an exciton in a parabolic quantum well as function of the elec- 
tron confinement strength Aj. The comparison is made for different approaches described 
in Ref. [12|. The parameters a 2 = m\jm,2 and x = ^2/Ai are fixed as indicated. 



The trial companions of the exact action ( |T4"D were combinations of oscillator 
trial actions for the center-of-mass and the ^-coordinate and three-dimensional 
(two-dimensional) Coulomb potentials for the three-dimensional (two-dimensional 
in-plane) relative coordinates. The corresponding results (see Ref. [0) can be found 
in the following figures and are denoted as quasi three-dimensional (Q3D) and quasi 
two-dimensional (Q2D or Q2Dalt) ansatz. In Fig. [| we neglect any phonon influ- 
ence to demonstrate the smooth interpolation of the limiting values 1-Rqo and AR^ 
of the binding energy (actually we plotted there the ground-state energy with the 
continuum edge being subtracted, that is, the quantity —Eb)- Fig. || shows results 
for the general case; we present data for the ground-state energy as well as the con- 
tinuum edge, which is the reference for the binding energy and has to be calculated 
separately. 
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FIG. 5. Ground-state energy of an exciton-phonon system in a parabolic quantum 
well as function of the confinement electron strength Ai. The remaining parameters 
r\ = yjRoo/hoj and £ = 1 — £oo/^o are fixed as indicated. In addition, an upper bound for 
the energy of the continuum edge is shown. 
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